Improved lower and upper bounds for entanglement of formation 
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We provide analytical lower and upper bounds for entanglement of formation for bipartite systems, 
which give a direct relation between the bounds of entanglement of formation and concurrence, and 
improve the previous results. Detailed examples are presented. 
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PACS numbers: 03.65.Ud, 03.67.Mn 

Quantum entanglement [l| is of special importance 
in quantum-information processing and is responsible 
for_many quantum tasks such as quantum teleportation 
, dense coding Q , swapping y| M , error correction 
and remote state preparation [H, [13] • The entan- 
glement of formation (EoF) [III E3 is a well-defined im- 
portant measure of entanglement for bipartite systems. 

Let Ha and Hb be m- and n-dimensional (to < n) 
vector spaces, respectively. A pure state \ip) € Ha <8> Hb 
has a Schmidt decomposition \ip) = Y^Li v^N)' wnere 
Hi > and 53i=i Mi = 1- The entanglement of formation 
is given by the entropy of the reduced density matrix 
p A =Tr B (\i ) ){^\), 

m 

E(\yj)) = S(p A ) = -£>log W s ff(/2), (1) 
»=i 

where log stands for the natural logarithm throughout 
the paper, pi (z = 1, 2, to) are the non-zero eigenvalues 
of pa and fx is the Schmidt vector (//i,//2, ...,p m ). For a 
bipartite mixed state p, the entanglement of formation is 
given by the minimum average marginal entropy of the 
ensemble decompositions of p, 



E(p) = min 



(2) 



for all possible ensemble realizations p — '^2 i Pi\ipi}('4 , i\y 
where pi > and ^ 4 pi = 1 . 

Another significant measure of quantum entanglement 
is the concurrence. The concurrence of a pure bipartite 
state \ip) is given by 



C(m = y/2[l - Tr(p\)} 



\ 



2(1 ( 3 ) 



It is extended to mixed states by the convex roof con- 
struction 



C{p)= min VftC(|^)), 



(4) 



for all possible ensemble realizations p = YliPi\''Pi)(' t Pi\- 

Considerable effort has been made to estimate the en- 
tanglement of formation and concurrence for bipartite 
quantum states, and their lower and upper bounds via 



analytical and numerical approaches. For the two-qubit 
case, EoF is a monotonically increasing function of the 
concurrence, and an analytical formula of concurrence 
has been derived [l3|. For the general high-dimensional 
case, due to the extremizations involved in the compu- 
tation, only a few analytic formulas have been obtained 
for isotropic states [l4| and Werner states [l5| for EoF, 
and for some special symmetric states for concur- 

rence. 

Instead of analytic formulas, some progress has been 
made toward the lower and upper bounds of EoF and 
concurrence for any m®n (to < n) mixed quantum state 
p. In [19|-[23j , explicit analytical lower and upper bounds 
of concurrence have been presented. In Ref. [19[ , a simple 
analytical lower bound of EoF has been derived. Recently 
new results related to the bounds of EoF have been fur- 
ther derived in [24|, • In this article, we give new lower 
and upper bounds of EoF based on the concurrence. De- 
tailed examples are presented, showing that our bounds 
improve the bounds in 0, [25( . 

Hj], the authors defined X(X) and Y(A) to de- 



In Ref. 



rive measurable lower and upper bounds of EoF. We give 
an improved definition of X(X) and Y(X) in this paper. 
For a given pure state = Yli=i \/Jk\ii)i the concur- 
rence of is given by c = y2(l — J2T I 1 !)- However, 
the pure states with the same value of concurrence c are 
not unique. Namely, different sets of the Schmidt coeffi- 
cients {pi} may give rise to the same concurrence c. The 
entanglement of formation H(p) defined in Eq.([T]) for a 
pure state depends on the Schmidt coefficients {pi}. We 
define the maximum and minimum of H(p) to be X(c) 
and Y(c) for fixed c, 



and 



X(c) = max ^ H(p) 



Y(c) = min ^ H(p) 



\ 



2(1 -5>?) 



\| 



2(1 -E^ 2 ) 



(5) 



(6) 



respectively, where the maximum and minimum are 
taken over all possible Schmidt coefficient distribu- 
tions {pi} such that the value of concurrence c = 
V^l-ErV?) ^ fixed. 
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Let e(c) be the largest monotonically increasing convex 
function that is bounded above by Y(c), and 77(c) be the 
smallest monotonically increasing concave function that 
is bounded below by X(c). 

Theorem 1 For any m ® n [m < n) quantum state p, 
the entanglement of formation E(p) satisfies 

e(C(p)) < E{p) < V (C( P )). (7) 

Proof. We assume that p = X^ftlV'iXV'il is the opti- 
mal decomposition of E(p). Therefore 

e(p) = E i pMi>i)) = E i PiHm 

> T,iPi £ ( c i) > £ (T,iPi c i) (8) 

> e(C(p)). 

We have used the definition of e(c) to obtain the first 
inequality. The second inequality is due to the convex 
property of e(c), and the last one is derived from the 
definition of concurrence. On the other hand, as p = 
SiftlV'iXV'il is the optimal decomposition of C(p), we 
have 

E{ P ) < Y.iPiE{m)=Y.iPiH{jk) 

< E 1 p^( c ») < ^(E ( k c (9) 



where we have used the definition of E(p) to obtain the 
first inequality. The second and the third inequalities are 
due to the definition of 77(c). □ 

Our analytical bounds ([7]) give an explicit relations 
between the EoF and the concurrence. In fact, if we 
denote c and c the analytical lower and upper bounds 
of concurrence, respectively, according to Theorem[T] we 
have the following corollary: 



Corollary 1 For any m (g) n (m < n) quantum state p, 
the entanglement of formation E{p) satisfies 



e{c) < E(p) < 77(c). (10) 



Here c and c could be any known analytical lower and 
upper bounds of concurrence. For example, from the 
known bounds of concurrence in EI [MSI, one may 
choose 



= n(C(o)), I 

c = max {0, ^ZgZ^dlpT. || _ 1); ^/ZgZ: [R(p) l] , ^[TtV) _ T r(p A )}, ^/2[Tr{p*) - Tr(p%)}} 



and c = mm{y/2[l-Tr(p 2 A )], y/2[l - Tr(p|)]|, where 

P Ta stands for the partial transpose with respect to the 
subsystem A, R{p) is the realigned matrix of p, \\ ■ \\ 
stands for the trace norm, and pa and pb are the reduced 
density matrices with respect to the subsystems A and 
B respectively. 

The maximal admissible H(p) and the minimal ad- 
missible H(p) in Eqs.© and ([6]) for a given c can be 
estimated following the approach in Ref. 24]. Let n\ be 
the number of entries such that m = a and let 77,2 be 
the number of entries such that pi = f3. The maximal 
admissible H(p) and the minimal admissible H(p) for a 
given c become, for fixed ni,n 2 ,ni + n 2 < m , one of 
maximizing or minimizing the function 



F„ 



), (11) 



where h(x) — — xlogx, 



+ \Jn\ - m{ni +n 2 )[l - n 2 (l - ^ )] 
ni{nx + n 2 ) 

and = (1 - nia nin2 )/n 2 . 

When to = 3, to find the expressions of upper and 
lower bounds in Eqs.((7]) and ([TO]) is to obtain the max- 
imization and minimization over the three functions 



Fu(c), F 12 (c) and F 21 (c). From Eq.(HI]), for m = 3, 
we have 

F n {c), 0<c<l, 

X ^ = < tp n 1 , < 2 (12) 

F 12 ( C ), 1< C< -y=, 



and 



Fu(c), 0<c<l, 

Y ^ = < J? M 1 ^ <r 2 (13) 
F 2i (c), 1< c < — . 

To determine e(c) and 77(c), we study the concavity and 
convexity of functions i'ii(c), i<i2(c) and F 2 i(c). Since 
F n > and F 21 < 0, where F^ is the second derivative 
of Fij, we have that e(c) is the curve consisting of Fn for 
c e (0,1] and the line connecting the points [l,i<2i(l)] 
and [-±=,F 21 (-^)} for c€ (1, -J=), that is, 

(F n (c), 0<c<l, 
e(c) = < %/31og3/2 . 2 

l^vr (c - 1)+l0g2 ' 1<c£ vT 

(14) 

Similarly, since F u > and F 12 > 0, we have that 77(c) is 
the curve connecting the points [0,0] and [-7=, ^2(77=)] 
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for c G (0, 1], and the line connecting the points [1, Fi2(l)] and [-^=, i*i2(^)] for c € (1, -7=], that is, 



log 2(c), 0<c<l, 
77(c) = <( 21og3/2 + log6-31og3 
3(^3-2) 



(\/3c-2) + log3, l<c<-p. 

v3 



(15) 



See Fig. 1. 
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FIG. 1. Upper and lower bounds of E(p) (dished lines and 
dotted lines ), and Fn,Fi2 and F21 (solid lines). 

Similarly, for any to, we can get the expressions of X(c) 
and Y(c), 



fFn(c), 0<c<l, 



F 12 (c), Kc< 



X(c) = I 



2(m-2) 
F 1{m -i)(c), r 1 <c< 



TO — 1 



2(m- 1) 



(16) 



and 



(F n (c), 0<c<l, 



F 2 i{c), Kc< 



y(c) = ^ 



V3' 



(m-l)l 



2(m- 2) 

TO — 1 



< C < 



2(m - 1) 



(17) 

The representations of e(c) and 77(c) can be also sim- 
ilarly calculated analytically in accordance with the fol- 
lowing principles (2 <t < m — 1): 



If F' u (c) > 0, c G (V^'V^' thCn 
+ F lt {JK); If K(c) < 0, c G 



2(t-l) 



0, ce(v /^ 

0, c G ( 
F n ( 



2(«-l) 



then 77(c) = F lt (c); If F t ' x (c) > 
then e(c) = F tl (c); If i£(c) < 
'^], then e(c) = [^(a/S) " 



The bounds given in Theorem 1 and Corollary 1 can 
be used to improve the bounds of EoF presented in 24 1 
and [25j]. In fact, the lower bound obtained in Ref. 25 j 
is better than the lower bound from Ref.[24|, while the 
upper bounds are the same. Our bounds are obtained 
from the improved bounding functions ([5]) and ^ . They 
are directly given by the concurrence. From the concur- 
rence, or the lower and upper bounds of the concurrence 
of a given mixed state, one can get analytical lower and 
upper bounds of EoF of the state. To see the tightness 
of inequalities (O and (fT0|). let us consider the following 
examples. 

Example. 1 Let us consider the well-known Werner 
states, which are a class of mixed states for d®d systems 
that are invariant under the transformations U ®U , for 
any unitary transformation U [H, H|| . The density ma- 
trix of the Werner states can be expressed as 



Pf 



1 



d 3 -d 



[(d - f)I + (df - 



(18) 



where T is the flip operator defined by T((j)(g)ip) = rp®<f>. 
Consider the case d — 3. We have (p/)a = {ps)b — 
I(|0)(0|+]l)(l| + |2)(2|) and 1 - Tr[(p f ) A ] = §. By 
Refs. [24] . [25l | . the upper bound of EoF is given by E(pf) < 
1.099. From Eq. (Tl5]) , we get the upper bound of pf, 



E( Pf ) < -/log 2, -l</<0, 



(19) 



It is obvious that -/log 2 < 1.099 for -1 < / < 0. Hence 
the upper bound (I15p is better than the upper bound in 
Refs.HlH. 

Example. 2 Consider the 3 ® 3 mixed state p = 



(1 — x)\ip) where the column vector \ip) 



2/ 

9 

(a, 0, 0, 0, ^, 0, 0, 0, Wa 2 + 2/3 with a G [0,1], * 

stands for vector transposition. For this state we have 
Tr(p 2 ) - Tr(p A ) = Tr{p 2 ) - Tr(p 2 B ) = 2 [9 - 26x + 
9a 4 (-2 + x)x + Ux 2 + 6a 2 (9 - 22x + llx 2 )]/9(2 + 3a 2 ) 2 
and 1 - Tr(p\) = 1 - Tr(p%) = [6 + 4a; - 18a 4 (-2 + 
x)x - 2x 2 + 12a 2 (3 - 2x + x 2 )}/3(2 + 3a 2 ) 2 . 



For fixed x = 0.1, one has 



c = v/2[Tr(p2) - Tr{p\)} 

2^6.53 + 41.46a 2 - 1.71a 4 
3(2 + 3a 2 ) : 



(20) 



c = v/2[l-Tr(^)] 

I J2(6.38 + 33.72a 2 + 3.42a 4 ) 
2 + 3a 2 ' 



(21) 



Substituting c and c into Eqs. (IT4l) and (|T5|) . we have the 
upper and lower bounds for E(p), see Fig. 2. 
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FIG. 2. Upper and lower bounds of E(p) when x=0.1. 
Dished lines are given by Eos. ((Til and (|15|) . solid lines are 
given by Ref. 0. 

Similarly for x = 0.001, one has c and c, and the upper 
and lower bounds for E(p), see Fig. 3. Form Figs. 2 and 
3, we see that the lower bound of EoF (fT4")l is better than 
the one from 12411. 





FIG. 3. Upper and lower bounds of E(p) when x=0.001. 
Dished lines are given by Eqs. lfH)) and (|15J) . solid lines are 
given by Ref. [241 ] . 

In fact, we can make the bounds better by choosing 
suitable c and c. For instance, to find a better lower 



bound, we may choose c = y^GI/ ^ 4 !! — !)■ Then for 
x = 0.1 and a £ [0.5,0.66], one has 

^ T { _ 2 [5 + 6.9a 2 - 0.9a 4 + 9.353a(2 + 3a 2 )] 



llP 11 ' 3(2 + 3a 2 ) 2 

When x = 0.001 and a € [0.57, 0.59], one has 



(22) 



T _ 2[5.99 + 8.978a 2 - 0.009a 4 + 10.382a(2 + 3a 2 )] 

l|P "~ ~ 3(2 + 3a 2 ) 2 ' 

(23) 

Substitute Eqs. ([22]l and (J23J) into (O, respectively, we 
get another lower bound of EoF for the state in example 
2, see Figs. 4 and 5. 
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FIG. 4. Lower bounds of E(p) when x=0.1. Dashed line is 
obtained by Ref. [25[, dotted line is obtained by Eq. lfTi")) and 
solid line is obtained by Ref. [24| . 
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FIG. 5. Lower bounds of E(p) when x=0.001. Dashed line 
is obtained by Ref. [23] . dotted line is obtained by Eq. (|14|l 
and solid line is obtained by Ref. [24l] . 



From Figs. 4 and 5, it is obvious that e(\/^(\\p" 

1)) > E\ > eo when x = 0.1 and a 6 [0.5,0.66], and 

ei = eW2[Tr{{?)-Tr{(? A )\) > £(^(11^11 - 1)) > s 

when x = 0.001 and a £ [0.57,0.59], where £o and e\ 
are obtained by Ref. 24] and Ref. 25] , respectively . The 
lower bounds are improved in the particular interval of 
a. 
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The density matrix in example 2 is close to being sepa- 
rable (pure) when x is close to 1 (0). To show better the 
advantage of our results, we now take x = 0.6. We have 



then Tr(p 2 ) - Tr(p 2 A ) = Tr(p 2 Y- Tr(p 2 B ) < 0, and hence 
the lower bounds from Refs.JU [25[ are e = £ i = 0. 
However, by choosing 



- 1 



2(6 + a 2 (9 - 21x) + 6-\/3a(2 + 3a 2 ) (x - 1) - lOx - 9a 4 x) 



3V3(2 + 3a 



2Y2 



(24) 



we have c > when a £ (0.205, 1). From Corollary 1 we 
get that E(p) > e(c) > 0, see Fig. 6 for the lower bound 

e(c), 




FIG. 6. Lower bounds of E(p) when x=0.6. 

In summary, we have presented analytic lower and up- 
per bounds of EoF for arbitrary bipartite mixed states. 
The bounds can be used to improve the previous results 
on bounds of EoF. Although the EoF is a monotonically 
increasing function of the concurrence only in the two- 
qubit case, it turns out that for higher dimensional cases, 
the bounds of EoF have a tight relation to the concur- 
rence or the bounds of concurrence. 
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